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Meson-meson interactions - from static to dynamic valence quarks* 
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A method for the extraction of an effective meson-meson potential from Green functions, which can be obtained 
from a lattice simulation, is presented. Simulations are carried out for compact QED and QCD in four dimensions 
using the quenched approximation and the hopping parameter expansion. In a further study, a heavy-light meson 
is considered employing a conjugate gradient algorithm for the light propagators. Due to the Pauli exclusion 
principle, the results for QED indicate the existence of a hard core, but for QCD there is strong attraction at 
small meson distances. 



1. INTRODUCTION 

QCD is believed to be the correct theory to 
describe the physics of strong interaction. There- 
fore it should be possible, in principle, to extract 
the potential between hadrons from this theory. 
Before we can proceed, we have to define what is 
meant by "potential" . Throughout this paper we 
will use two different definitions: 

1. The potential is the function V{r) that has 
to be inserted into a Schrodinger equation 
so that quantum mechanics and quantum 
field theory result in the same predictions. 

2. The potential at finite temperature is ap- 
proximated by the thermodynamic free en- 
ergy F{r). 

For the first definition we will carry out a simu- 
lation at zero temperature and demand that this 
simulation and the solution of a corresponding 
Schrodinger equation give the same Green func- 
tions for the two-meson system. For the sec- 
ond definition we will follow the guidelines of the 
finite-temperature concept [1], generalized to dy- 
namic valence quarks. 

2. POTENTIAL AT T = 

We wish to compare the Green functions from 
quantum mechanics with those from field theory. 
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In a euclidean formulation of quantum mechanics 
with imaginary time r the 4-point Green function 
representing the two-meson system can be written 
as 

G^'^ = Y,{yuy2m^i,^2)e-'^-'' . (1) 

n 

This describes the propagation of two mesons 
from xi,X2 at T = to t/i,t/2 at r = Tf . E„ 
are the energy eigenvalues and \n) the associated 
eigenstates. For a given potential we can solve 
the Schrodinger equation to obtain E„ and \n) 
and calculate the sum in (1). In order to invert 
this process and calculate the potential from a 
given Green function we have to make some ap- 
proximations and accept some restrictions. The 
first step is to set the final positions y equal to 
the initial positions x. It leads to 

G^^\d,Tf) = J2\{^i,^2\n)\'e-^'^^', (2) 

n 

with a particle separation d = X2 — xi. In the 
next step we neglect all terms of the sum except 
the largest contribution, which is assumed at, say, 
n = m for a given d: 

G'^^\d,Tj)^\{x,,X2\m){'e-^-''. (3) 

The question emerges how Em and [xi, X2\m) are 
related to the potential V{x2 — xi). We try to es- 
timate the magnitude of the terms in the sum. 
Each term is a product of the squared modulus 
of the relative wave function and e~^"'^f . We can 
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Figure 1. Estimate for the terms contributing 
to the sum in the spectral decomposition of the 
Green function (2). 



neglect the term if either the wave function or the 
exponential is small. Because the wave function 
decays exponentially "outside" the potential, it 
can be neglected at a given d when E„ < V{d). 
On the other hand, e'^^"^' goes to zero if EnTj 
becomes large and for fixed tj the larger E^s 
need not be considered. The maximum contribu- 
tion for a given d and tj can be expected near the 
point En ~ y{d)- In this sense the eigenvalue E^ 
is an approximation for the potential (see also fig. 

1): 

G(4)(rf,r;) « c(rf)e-^('')^^ (4) 

with the positive factor c[d) = a;2|"^)P- 
From this relation we obtain the potential as 

« (inc(rf) -lnG(4)(rf,r;)) . (5) 

3. POTENTIAL AT T 7^ 

In the case of finite temperature we take the 
thermodynamic free energy [1]: 

F = -T\nZ, (6) 



with the temperature T = and the partition 
function 

Z = Y.{s\e-^"\s), (7) 

s 

where s contains the external sources and the vac- 
uum states. The two-meson wave function in pe- 
riodic time 'iI>[xi,X2,t) evolves from 'il>[xi,X2,Q) 
according to the Green function: 

i){xi,X2, t) = 

j dyidy2G^'^^xi,X2,T;yi,y2,0)ip(yi,y2,0). (8) 

Proceeding as in [1], we obtain the partition func- 
tion 

Z = {G^''Hxi,x2,^;xi,x2,0)). (9) 

If we assume in F{d) = V{d) — TS{d) that the 
entropy depends only smoothly on the separation 
of the external sources, we can approximate the 
potential by the free energy. 

4. GREEN FUNCTIONS FROM QCD 

In the preceding sections we found a corre- 
spondence between the 4-point Green function for 
mesons and the potential in quantum mechan- 
ics. These quantum mechanical 4-point functions 
represent 8-point functions in QCD, because each 
meson consists of a quark and an antiquark. The 
additional colour degrees of freedom have to be 
saturated in a gauge invariant way. At zero tem- 
perature this is done as 

Gf]^bccdd {yi,yi,y2,y2,T;xi,xi,x2,x2,o), (10) 

meaning that quark and antiquark are contracted 
at the endpoints of the 8-point function. At finite 
temperature the contraction is made over the pe- 
riodic time: 

G(^^^^-^°^{x,,x2,^;x,,x2,0) = 

^Edcdl^i' *2, a;2, «i, «i, «2, «2, 0). (11) 

The 8-point functions G'-®-' can be split into 
2-point functions G'-^-' (quark propagators) by 
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Figure 2. Meson-meson potential as a function of relative distance d. For U{1) a short-ranged hard core 
is observed due to antisymmetrisation of the propagator. In the case of SU{3) an attractive interaction 
is obtained since the nonabelian colour structure does not lead to a vanishing propagator. 



fermionic integration in the path integral. This 
gives the same result as the application of Wick's 
theorem. Thus we can write [2]: 

G'^^\xi,Xi,X2, X'j, T; Xi, Xi,X2, X2, 0) = 

A + B - C - D, (12) 

with the term 

A = G^'^^xi,t;xi,0)G''^'^^xi,t;xi,0) 

G^^\x2,t; X2, 0)Gt(2)(a;2, r; X2, 0) (13) 

representing the direct propagation of all quarks 



and antiquarks from {xi,0) to {xi,T), B repre- 
senting the exchange of quark-antiquark pairs, 
and G and D reflecting the exchange of one quark 
or antiquark. The 2-point functions can be ob- 
tained from a lattice computation. 

5. RESULTS 

The numerical simulation was carried out on 
a lattice with spatial extension x Ny x = 
8x8x16 with periodic boundary conditions. The 
temporal extension was chosen to be Nr = 16 for 
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zero temperature and Nr = 4 for finite temper- 
ature. The gauge fields were generated in the 
quenched approximation and 200 configurations 
separated by 50 sweeps each were analyzed. For 
the SU (3) simulation the inverse coupling was set 
to /3 = 5.6 and for U{1) to /3 = 0.9. The propaga- 
tors were calculated for Kogut-Susskind fermions 
using the hopping parameter expansion. For con- 
vergence in this expansion, the quark mass was 
set to niqa = 8. 

Fig. 2 shows the results for the extracted po- 
tentials. We assume a weak rf-dependence of 
c{d) « const. The cluster value was subtracted 
to give V{d co) = 0. The potentials for U{1) 
exhibit a hard core both for zero and finite tem- 
perature (left hand side of fig. 2) which can be un- 
derstood analytically. At c? = quark exchange 
and direct terms cannot be distinguished and the 
terms A through D in (12) are equal. The 8- 
point function is thus zero and the potential in- 
finite. For SU{3) there is strong attraction at 
d = (right hand side of fig. 2). This happens be- 
cause the particle exchange results in a different 
saturation of the colour degrees of freedom and 
the terms A through D cannot cancel each other. 
For meson distances c? > 1 all potentials indicate 
negligible interaction compared to the statistical 
error bars. This narrow range can be expected 
because of the small meson size which is inversely 
proportional to the quark mass chosen very large 
for this simulation. 



6. HEAVY-LIGHT SYSTEM 

In a second approach we inverted the fermionic 
matrix by a conjugate gradient algorithm using 
random source techniques. One of the quarks in 
the meson was chosen static and the light quark 
mass was set to niqa = 0.1 [3]. An 8 x 8 x 8 x 16 
lattice with 100 SU{3) configurations at /3 = 5.6 
separated by 200 sweeps was considered. The 
potential is extracted using (4), however the d- 
dependence of c{d) is not neglected now. An ex- 
ponential fit to the time behaviour of the Green 
function is made, with c{d) and V{d) as param- 
eters. The resulting potential is drawn in fig. 3 
and it turns out to be attractive again. 
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Figure 3. Meson-meson potential for heavy-light 
systems in SU{3) at T = 0. Attraction is ob- 
tained as for the SU{3) results in fig. 2. 



7. CONCLUSION 

We found a qualitative difference between the 
meson-meson potentials in compact QED and in 
QCD for small meson distances. In QED there is 
a hard core which can be understood analytically 
as the cancellation of the terms in a Wick expan- 
sion [2]. In QCD this cancellation is not enforced 
due to the sum over the colour degrees of freedom 
and no hard core is seen. However, all simulations 
were carried out in simplified models: with an un- 
realistically high quark mass in one case, with one 
quark of the meson static in the other. Our next 
aim is to investigate the full dynamic system with 
more realistic quark masses, using the conjugate 
gradient algorithm. 
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